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The negatively charged nitrogen- vacancy (NV~) color center in diamond has generated much 
interest for use in quantum technology. Despite the progress made in developing their applications, 
several questions about the basic properties of NV~ centers remain unresolved. In particular, 
knowledge of the phonon modes of the ^ Ai electronic state is key for understanding the optical 
pumping process. Using pump-probe spectroscopy, we measured the phonon sideband of the ^E 
^ Ai_ electronic transition in the NV~ center. From this we calculated the ^E — > ^ Ai_ one-phonon 
absorption spectrum and found it to differ from that of the '^E — >■ A2 transition. 

PACS numbers: 78.40.Ha, 63.20.kp, 63.20.Pw, 61.72.jn 

Keywords: Nitrogen-vacancy centers, phonon-defect interactions, excited-state spectroscopy 



The nitrogen-vacancy (NV) center in diamond 
(Fig. [1^) is a color center consisting of a substitutional 
nitrogen atom in the diamond crystal lattice adjacent to 
a missing carbon atom (a vacancy). NV centers have 
C31, point-group symmetry and have discrete electronic 
energy states between the diamond valence and conduc- 
tion bands. The negatively charged NV~ center can 
be optically spin-polarized and read out, and it has a 
long ground-state transverse spin relaxation time at room 
temperature [TJ |2] • These properties make NV~ centers 
useful in a variety of applications including electric and 
magnetic field sensing |3H5] , rotation sensing [7H3] , quan- 
tum computing [10' '1 II , quantum cryptography f l2[ I13j , 
and sub-diffraction-limited imaging [14h16j. Despite the 
progress made on developing these applications, the com- 
plete NV~ energy level structure and vibronic structure 
are unknown. 

Figure [ijD shows a simplified NV~ energy-level dia- 
gram as confirmed by experiment. The triplet-triplet 
(^^2 O ^E) and singlet-singlet {^E o ^Ai) energy dif- 
ferences are known to be 1.945 eV (637 nm) and 1.190 
eV (1042 nm), respectively [T7H^ . However, neither 
where these energy states lie with respect to the di- 
amond valence and conduction bands nor the triplet- 
singlet ('^^2 O and ^i? O ^E) energy differences 
are known directly [21] [22] . Theoretical calculations pre- 
dict the existence of additional energy states {}E' and 
^A'^), but disagree on their energies (see Refs. [^5H?7] 
and references therein). Prior experiments and ah initio 
calculations studied the phonon sidebands (PSBs) for the 
3^2 and ^E transitions [T7|l28H3T]. The 

^E -> ^Ai and ^Ai -> ^E PSBs have not been studied 
theoretically, and only the ^Ai — >■ ^E transition had been 



measured prior to this work [T5H20j . 

A more complete experimental picture of NV^ proper- 
ties can provide insight for applications and validate the- 
oretical models of NV^ center attributes. We attempt 
to fill the gaps in the knowledge of NV~ properties by 
measuring the ^E — ^ ^Ai PSB and searching for pre- 
viously unobserved transitions. Finding the ^E ^E' 
zero-phonon line (ZPL) would resolve the disagreement 
on the predicted energy. The ^E — > ^Ai PSB yields 
information about the ^Ai phonon modes, which are also 
of interest. The spin-orbit interaction mixes the ^E and 
^Ai states, resulting in triplet-singlet intersystem cross- 
ing (ISC). This enables non-radiative decay from the 
nominally '^E state to the nominally ^Ai state. The 
ISC rate is comparable to the '^E ^A2 spontaneous 
decay rate |32L I33| and is an important factor in the op- 
tical pumping process to the ^A2 rus ~ state. Mea- 
suring the ^Ai phonon modes could provide insight on 
the mechanism used for NV" spin polarization and read- 
out. Furthermore, the accepted group theoretical model 
of NV" predicts '^^2 and ^Ai to have the same elec- 
tronic configuration, meaning they should have the same 
phonon modes. A comparison between the ^A2 ^E 
and ^E — > ^Ai PSBs should be sensitive to differences 
between the ^^2 and ^Ai configurations. 

In our experiment, we populated the metastable ^E 
state using pump-laser light and measured transmission 
of probe-laser light through a diamond sample contain- 
ing an ensemble of NV~ centers (Fig. [ija and Fig. |2|. 
We determined the probe transmission through the di- 
amond with and without NV~ centers in the ^E state. 
A 300 mW 532 nm frequency-doubled Nd:YV04 pump 
laser beam and a 5 mW supercontinuum probe laser 
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FIG. 1. (a) The diamond lattice structure, containing an NV 
center, (b) The NV~ energy-level diagram and our pump- 
probe spectroscopy scheme. The states are labeled by their 
Czv representations and electron spin multiplicities. Solid ar- 
rows are optical and microwave transitions, and dashed arrows 
are non-radiative transitions. The label "ISC" indicates inter- 
system crossing, which occurs primarily for the rus — ±1 
states and is responsible for optical pumping, (c) A config- 
uration coordinate diagram for Ai phonon modes showing 
the harmonic nuclear potential wells and phonon energy lev- 
els. The configuration for each electronic state is denoted in 
parentheses, and Qai is the normalized nuclear coordinate. 
The and states should have afe'^ configuration, 
should have flie^, and ^E should have a mixture of the two 



beam (wavelength range 450-1800 nm) were combined 
on a dichroic beamsplitter and focused with a 40x mi- 
croscope objective onto a cryogenically cooled diamond 
sample. The transmitted light was collimated and de- 
tected with a spectrometer with ~1 nm resolution. A 
chopper wheel modulated the pump light and a com- 
puter collected a transmission spectrum each time the 
pump light was blocked and unblocked. Absorption from 
appeared as a difference between the "pump blocked" 
and "pump unblocked" supercontinuum transmitted in- 
tensities. In another experiment, we used 912 nm and 
1042 nm continuous-wave (cw) lasers as probe sources 
and replaced the spectrometer with a photodiode 

Figure [Sji shows the '^E -J> '^Ai ZPL and PSB super- 
continuum absorption spectrum taken at 10 K with the 
sample "B8", a synthetic type lb high-pressure, high- 
temperature (HPHT) diamond with ^10 ppm NV~ con- 
centration. The PSB includes narrow absorption lines at 
811 and 912 nm and broad absorption features at 872, 
922, 931, and 983 nm. We believe the 811 and 912 nm 
lines are due to a 169.28(4) meV phonon mode and that 
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FIG. 2. The experimental apparatus. The data acquisition 
device (DAQ) monitors the chopper wheel state and trig- 
gers a spectrum acquisition when the pump is blocked or un- 
blocked. The computer collects "pump blocked" and "pump 
unblocked" transmission spectra. 



the other lines are due to a distribution of phonon modes. 
Figure shows the ^E — >■ "^^2 fluorescence spectrum 
taken at 4 K with a similar diamond. This PSB has a 
broader energy range, and has features at 686, 692, and 
696 nm. Using these measured spectra and the tech- 
niques outlined in Refs. [551 IMl 131], we calculated the 
^E ^Ai and ^E ^A2 one-phonon spectra (Fig. |4|, 
which are the rates at which these transitions create one 
phonon of a given energy. We expect these one-phonon 
spectra to be similar, since both come from E ^ A tran- 
sitions with the same final-state electronic configuration 
(Fig. [ij;) . The one-phonon spectra show rough similarity, 
and the differences between them are because of higher- 
order corrections to the ^Ai and '^^2 electronic states. 

We observed the above ^E — >■ ^Ai PSB features in sev- 
eral diamond samples, and the absorption was greater in 
samples with higher NV~ concentration. The ^E — >■ ^Ai 
absorption should increase with pump power and satu- 
rate when the pumping rate becomes comparable to the 
^E decay rate. The absorption at room temperature in- 
creased linearly with pump power (up to 60 mW focused 
to a beam waist smaller than 5 \xm), indicating that the 
^E population was not saturated. However, the absorp- 
tion at 10 K saturated at ~15 mW. This saturation is 
likely due to the enhanced ^E lifetime at cold tempera- 
ture [18] . Introducing a static transverse magnetic field 
to the samples improved the absorption contrast by a few 
percent. This is because the Zeeman interaction mixes 
the triplet spin sublevels, which spoils the optical pump- 
ing to rus = and increases the ^E population. We did 
not detect ei^E ^E' ZPL in the 480-1100 nm range of 
the supercontinuum transmission spectrum, which means 
this transition lies outside of this range or was too weak 
to detect. This wavelength span was limited by the spec- 
trometer. 

Using a rate equation calculation based on the NV~ 
excitation and decay rates at room temperature [33], we 
estimate the room temperature ^E — >■ ^Ai ZPL cross 
section to be roughly 4 x 10~^^ m^ [31], which is consis- 
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FIG. 3. (a) The supercontiuum absorption spectrum collected 
at 10 K for diamond sample B8 using 35 mW of pump-laser 
light. PSB fluorescence from — >■ '^A2 is present for wave- 
lengths shorter than 840 nm and has been subtracted out. 
The horizontal lines indicate the expected PSB absorption 
energies for 71 and 169 meV phonons. (b) The fluorescence 
spectrum of a similar diamond collected at 4 K. The horizon- 
tal lines indicate the expected PSB absorption energies for 
64 meV phonons. Although the 686, 692, and 696 nm fea- 
tures are often ignored, they are vital to our comparison of 
the — 5> ^Ai and ^E — > PSBs, as they give rise to peaks 
(3)-(5) in Fig.|4] From these plots, we calculated Huang-Rhys 
parameters of 0.9 (a) and 3.49 (b). 



tent with previous work |36| . The accuracy of this cross 
section estimate is primarily limited by uncertainty in 
the NV^ concentration; varying the NV^ concentration 
from 5 to 20 ppm in our model yields estimated cross 
sections ranging from 3.4 to 5.4 x 10~^^ m^ (compared 
to 4.0 X 10-22 m^ with 10 ppm NV"). 

We varied the temperature of sample B8 from 10 to 
300 K and recorded the absorption-feature contrasts, 
linewidths, and integrated areas |34j- The features be- 
come weaker and broader with increasing temperature, 
and their integrated areas decrease. This decrease in 
area is consistent with the lifetime decrease observed 
in Ref. HSI. 



-e 2.5 



2.0 



c 1.5 



1.0 



.5 0.5 

T 

























(2) (3) (4) (5) ■ 




/e^Mi \ 







m 

2.0 J-^ 
o 

1.5 ^ 



1.0 



0.5 o 



- 



50 100 
Phonon energy (meV) 



150 



FIG. 4. The one-phonon spectra for the ^E — 5> ^Ai and ^E ^ 
transitions, extracted from Fig.js] The above spectra are 
normalized to have equal areas, and the ^A2 curve is vertically 
offset for clarity. In each spectrum we see five peaks, labeled 
(l)-(5), though the ^Ai peaks are shifted to higher energies 
(see Tab.|l]|. 



Peak # 


state 


^Ai state 


(1) 


64 meV 


71 meV 


(2) 


122 meV 


125 meV 


(3) 


138 meV 


141 meV 


(4) 


153 meV 


156 meV 


(5) 


163 meV 


169 meV 



TABLE I. The energies of the one-phonon peaks shown in 
Fig. [4] When comparing the energies of the ^^2 and ^Ai 
phonon modes, we see a systematic shift to higher energy of 
a few meV. 



Using cw probe lasers and a similar diamond sample 
"S2" (16 ppm NV" concentration), we measured the cen- 
ter wavelengths of the 912 and 1042 nm absorption lines 
at 40 K to be 912.19(2) nm and 1041.96(2) nm. Our 
ZPL center wavelength is consistent with previous mea- 
surements [T5J \W\ . At low temperatures, the 1042 and 
912 nm features have narrow widths (currently limited 
by the spectrometer resolution). These narrow widths 
imply that the vibrational mode associated with the 912 
nm feature is sharp. By measuring 912 nm absorption as 
a function of light polarization angle, we found that the 
912 nm absorption has the same polarization selection 
rules as the ^Ai ZPL [HIHl]. 

Comparing the ^E ^Ai absorption PSB depicted in 
Fig. [3^ with previous observations of the ^Ai ^E fluo- 
rescence PSB [inilSn], it is evident that these PSBs differ 
significantly. This difference is due to the anharmonicity 
of the ^E vibronic levels induced by the dynamic Jahn- 
Teller effect, which is not present in ^Ai |20]. In the 
low-temperature limit, the PSB features of A ^ E elec- 
tronic transitions will exhibit anharmonicity, while the 
PSB features oi E —i' A transitions will be harmonic [37] . 
Consequently, it is appropriate to compare the ^E ^ ^Ai 
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absorption PSB with the — > fluorescence PSB. 
Furthermore, ^Ai and have the same electronic con- 
figuration (of e^), and hence should have the same nuclear 
equilibrium positions and phonon modes. Since their ini- 
tial states are difi^erent, the ^Ai and ^E — >• 
transitions may couple to a different number of phonon 
modes and have different Huang-Rhys parameters, but 
the ^Ai and '^A2 one-phonon spectra should be the same. 

As mentioned above, we calculated the one-phonon 
spectra from the PSBs shown in Fig. [3j The n-phonon 
spectrum is the convolution of the (n-l)-phonon and one- 
phonon spectra, and the sum of all n-phonon spectra 
generates the transition PSB. The one-phonon spectra 
are also related to the ^Ai and '^^2 phonon density of 
states (DOS). As seen in Fig. [4j we found similarities 
between the one-phonon spectra; both spectra have one 
large feature and four small features. However, all of 
the ^E — >• ^Ai features are displaced to higher energies 
(Tab.|l]). 

Introducing a point defect into a lattice alters the vi- 
brational motion of the defect and its neighbors from 
what it would have been with ordinary atoms in the 
lattice. This is because the parameters that determine 
the frequencies of the vibrational motion for these atoms 
(the masses and effective spring constants) are modified. 
When the frequencies of the local oscillations of the de- 
fect lie within the spectrum of allowed vibrational modes 
of the remaining crystal, the local modes hybridize with 
the lattice modes and are called "quasilocal" (quasilocal 
because the nuclear oscillation amplitudes fall off slowly 
with increasing distance from the defect) [301 [35]. The 
^^71 meV phonon modes we observed appear to be from a 
quasilocal mode of NV^ in the ^Ai state. The diamond 
lattice phonon DOS is appreciable at 71 meV [551 139j. 
and since the NV~ 71 meV mode couples strongly to the 
diamond lattice modes, the peaks of the 71 meV mode 
are consequently broadened. 

In contrast to the quasilocal mode case, when the fre- 
quency of the local oscillations of a defect lies outside the 
lattice phonon DOS, this yields a "localized" mode. In 
this instance, the oscillations of the defect couple poorly 
to the oscillations of the rest of the crystal, the vibra- 
tional motion is confined to the region of the defect, and 
the local phonon mode energy is unbroadened. This is the 
case for the 169 meV mode. The diamond lattice phonon 
DOS has an upper limit of 168 meV [SHO]. The NV" 
169 meV mode falls outside the diamond lattice phonon 
spectrum and couples poorly to the lattice modes, conse- 
quently making the peaks of the 169 meV mode in Fig. [3^ 
sharp. 

The existence of a 169 meV local phonon mode and 
the differences between the ^E ^Ai and '^E — > •^^2 
one-phonon spectra are surprising for several reasons. 
Ab initio calculations for the NV~ triplet-state vibra- 
tions do not predict the existence of high-energy local 
phonon modes [IHl ISO], and the '^E ^ ^Ai PSB is the 



only NV~ PSB to contain such a feature. Due to the 
discrepancy in one-phonon spectra, we conclude that the 
^Ai and ^^2 levels have corrections to their lowest-order 
electronic states that affect the phonon modes of these 
states differently. Since the features in the one-phonon 
spectrum are shifted to higher energies, we can conclude 
that the nearby atoms are more tightly bonded in the 
^Ai state than in the '^A2 state. 

In summary, we measured the ^E — >■ ^Ai absorption 
spectrum of the NV~ center using pump-probe spec- 
troscopy. In the ^E — > ^Ai PSB and one-phonon absorp- 
tion spectrum we found several phonon modes, one of 
which lies outside the diamond lattice phonon DOS. The 
^E — ^Ai and "^E — ^A2 one-phonon spectra show gen- 
eral similarities, but the ^Ai phonon modes are shifted 
to higher energies, which is due to corrections to the ^Ai 
and ^A2 orbital configurations. We also searched for the 
^E ~^^E' ZPL for energies up to 2.0 eV, but we did not 
detect it. Since the ^E — >■ ^Ai and ^E — > ^E' transi- 
tions should have similar oscillator strengths, this sug- 
gests that the — > ^E' ZPL energy is greater than 2.0 
eV. Follow-up experiments will extend the search for the 
^E ^E' ZPL to higher energies with improved sensi- 
tivity. 
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K. Jensen, N. Leefcr, P. Hemmer, N.B. Manson, and D. Budker 
(Dated: January 25, 2013) 

INTRODUCTION 

In this supplement we present additional content to the results described in the main 
article. We include additional experimental technical details and findings, a description of 
how we calculated the one-phonon spectra from the — > ^Ai and — > absorption 
and emission PSBs, and an analysis of the resulting and ^Ai vibronic bands. 

EXPERIMENTAL DETAILS 

We used an Olympus LUCPlanFL N 40x PH2 microscope objective to focus the pump 
and probe beams onto the diamond samples (to a beam waist smaller than 5 (xm). This 
objective lens is achromatic, which ensures that the pumped NV" and probe spatial regions 
overlap. We used a Janis ST-500 liquid-hehum flow cryostat for cooling the diamond sample. 
Our pump-lascr source was a Coherent Verdi- V6, and our cw probe-laser sources included 
2 mW of 912 nm light from a Coherent CR 899 Ti: Sapphire laser, 30 mW of 912 nm light 
from a diode laser (1 nm linewidth), and 1.5 mW of 1042 nm light from an external-cavity 
diode laser (ECDL). We used a Fianium SC450-2 supercontinuum laser as our broadband 
probe and a Ocean Optics USB2000-I-VIS-NIR (~1 nm resolution, optimized for infrared 
sensing) for detecting transmitted supercontinuum light. 

We used the ^E — )■ ^A2 fluorescence spectrum of a similar diamond sample at 4 K for 
comparison with our ^E — ?■ ^Ai spectrum. This sample was also illuminated with 532 
nm pump light, and the emitted fluorescence was dispersed with a monochromator (0.1 nm 
resolution) and detected with a cooled-Ge detector (calibrated with the blackbody spectrum 
from a 3100 K tungsten bulb). Although we compared spectra from different diamond 
samples, the spectra from different high-concentration samples are in general consistent. 
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Ml ZPL CROSS SECTION ESTIMATE 

Using decay-rate and pump-laser absorption cross section parameters determined in [1, 
2], we constructed a rate equation model to estimate the fraction of NV~ centers in the 
metastable state throughout the diamond sample, from which we estimated the — > ^Ai 
ZPL cross section at room temperature. This model takes into account the pump beam 
divergence and absorption in the diamond. Using this calculation of the metastable NV~ 
center density and the experimentally determined 1042 nm transmission in sample B8 at 
various pump powers, we determined the — > ^Ai ZPL cross section. For simplicity, 
we approximated the probe beam to be a straight line through the pump beam axis. We 
estimated a — > ^Ai ZPL cross section of 4 x 10~^^ m^. This estimate is consistent with 
that of Ref. [3]. As mentioned in the main text, the uncertainty in NV~ center concentration 
dominates the error on this cross section estimate. Uncertainty in the NV^ center excited- 
state decay rates [1], the pump-laser absorption cross section [2], the pump beam waist, 
and the distance between the pump beam focus and the diamond surface contribute an 
additional uncertainty of about 0.9 x 10~^^ m^. 

TEMPERATURE DEPENDENCE OF THE Mi ZPL AND PSB 

We measured ^E — )■ ^Ai supercontinuum absorption spectra while varying the temper- 
ature of sample B8 from 10 to 300 K. Figure SI shows the absorption, linewidths, and 
integrated areas for the 1042 nm ZPL and the 811, 912, and 983 nm PSB features. These 
features became weaker and broader with increasing temperature, and the integrated areas 
decreased. The integrated areas should be independent of temperature [4]. However, the ^E 
state has shorter lifetime at higher temperature due to the enhanced electron-phonon decay 
rate to [5]. We believe the decrease in integrated area is because of the consequent 
reduction in ^E population at higher temperature. 

912 NM POLARIZATION SELECTION RULES 

We investigated the light-polarization selection rules for 912 nm absorption and compared 
them to those hsted in Table SI. An E ^ Ai transition is dipole-allowed for (x, y)-polarized 
light, while an E ^ E transition is also dipole-allowed for 2;-polarized light. A difference 
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FIG. SI. The percent absorption, linewidth, and integrated area oi^E ^ ^Ai absorption features in 
sample B8 as a function of temperature. The 811 and 912 nm lines vanished in the super continuum 
absorption spectra above 100 and 200 K, respectively. The spectrometer resolution contributes to 
the apparent linewidths plotted above; the gray region (3 meV) indicates where the spectrometer 
instrument broadening significantly contributes to the measured linewidths. We believe the 50 
K and 90 K measurements to be outliers due to thermal expansion in the cryostat during the 
measurement. The above error bars are one-sigma statistical errors extracted from the parameter 
fits of the absorption spectra. We estimate a 1 meV systematic uncertainty on the above linewidths. 



between the 912 nm selection rules and the expected E Ai ZPL selection rules could 
indicate that the 912 nm line is an £^ ^ £^ transition or that the ^Ai selection rules 

are not strictly obeyed in PSB transitions [5] . We determined the polarization dependence of 
912 nm absorption in a room-temperature optically-detected magnetic resonance (ODMR) 
experiment with diamond sample S2. We singled out the [111] -oriented NV" centers with 
an axial 15 G static magnetic field, exposed the sample to microwaves from a nearby wire, 
and measured the diode-laser absorption as a function of microwave frequency. Microwaves 
resonant with = ^ = ±1 transitions spoil the optical spin polarization, 

increase ^E population, and enhance probe absorption. By measuring the = — > 
nis — ±1 absorption in [lll]-oriented centers as a function of polarization angle for probe- 
light wavevector k parallel and perpendicular to the [111] z-axis, we found that the 912 
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nm transition is {x, ?/)-allowed and 2;-forbidden (Fig. S2) . These selection rules indicate 
an — 7- y4i transition and are consistent with the — ?► ^Ai ZPL selection rules [5]. We 
performed this experiment at 40 K with the 912 nm Ti: Sapphire laser and obtained consistent 
results for k \\ z, but we were unable to test k z because of mechanical constraints. 



Transition 


Dipole- allowed 


Ai ^ Ai 


z 


Ai ^ A2 




Ai^ E 


x,y 


A2 A2 


z 


A2<^ E 


x,y 


E ^ E 


x,y,z 




TABLE SI. Photon polarizations for dipole-allowed transitions between C^y electronic states [4]. 
The notation "x, y" implies that any polarization in the x-y plane has the same transition ampli- 
tude. The drawing on the right indicates the choice of axes. 




FIG. S2. (Left) A sample ODMR spectrum (k _L z) taken at room temperature, with absorption 
peaks for [lll]-oriented centers labeled. The middle two peaks are due to the other three orien- 
tations. (Right) 912 nm absorption by [lll]-oriented centers at different light polarization angles. 
The II z case is only sensitive to x,y polarization while the k J- z case is also sensitive to z 
polarization. The constant nonzero k \\ z ODMR peak height and sinusoidal variation to in the 
k J- z data indicate the transition is {x, y)-allowed and z-forbidden. 
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VIBRONIC ANALYSIS 

The vibronic structure predicted by the current NV~ model 

By adopting the adiabatic and harmonic approximations, the nuclear vibrational motion 
associated with a given electronic state is approximately governed by the harmonic potential 
formed by the dependence of the state's electronic energy on the nuclear coordinates. The 
nuclear vibrational potential associated with the n*'* electronic state is [6] 



En{Q) = En{0) + J2 (^n,a,j,kQa,j,k + -j{^a,j,k + K,a,j,k)Q a,j,k 
a,j,k 

1 



En{0) - SEn{Q) + 2 E i'^l,j,k + bn,a,j,k){Qa,j,k " SQn,a,j,ky 
a,j,k 



(1) 

(2) 



where 



En{0) = 




(^n,a,j,k — 




^a,j,k ~ 




bn,a,j,k — 






E 




a,j,k 



9Qa,j,k 
'^Qa,j,k 



1$. 



I$n 



dQi 



\^n)-(^L,k 



a. 



n,a,j,k 



0'n,a,j,k 



2(<^a,,-,fe + bn,a,j,k) 



(3) 



In the above expressions, He is the electronic Hamiltonian, Qa,j,k is the normal nuclear 
displacement coordinate (with respect to the nuclear equilibrium coordinates of the ground 
n = electronic state) of the a*'* ground electronic state eigenmode with symmetry [j, k) = 
{Ai, A2, (E.x), {E,y)} and vibrational energy huja,j,k- Note that quadratic terms that 
couple the eigenmodes in excited electronic states have been ignored and that for the ground 
electronic state ao,a,j,k = is a condition of nuclear equilibrium. It is clear that linear 
interactions displace the nuclear equilibrium coordinates by 5Qn,a,j,k and that the quadratic 
interactions shift the vibrational energies by Sujn,a,j,k ~ bn,a,j,k/'^- 

The vibrational state |xn,i/) associated with the n^^ electronic state with vibronic energy 
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En,iJ follows trivially from the vibrational equation [T{Q) + En{Q)] \Xn,u) = En,u \Xn,u), where 
T[Q) is the nuclear kinetic energy. Importantly, the vibronic coupling of the adiabatic 
vibronic states \^n,Xn,u) by Jahn- Teller interactions between degenerate electronic states 
has not been considered thus far. The Jahn- Teller interactions occur between degenerate 
electronic states and degenerate phonon modes and result in anharmonicities in the vibronic 
structures of the degenerate electronic states. For the case of the NV~ center, only E- 
symmetric electronic states and phonon modes exhibit degeneracy. As discussed in the 
article, the vibronic couplings and anharmonicites induced by the Jahn- Teller effect are not 
relevant to the PSBs analyzed in this work. 




FIG. S3, (color online) Configuration coordinate diagrams for Ai (Qai) and E {Qe) phonon modes 
depicting the harmonic nuclear potential wells and phonon energy levels. The optical emission and 
IR absorption transitions (solid arrows) are depicted according to the Franck-Condon principle. 
The current model predicts that the wells of the ^E and are displaced from the equilibrium 
coordinates (as indicated by dashed lines) differently, but the well of ^Ai is not displaced. The 
displacements of the degenerate electronic levels oi^E and ^E are equal and opposite (displacements 
do not imply static Jahn- Teller effects). The current model also predicts that the phonon energies 
of ^E and ^E', but not ^^i, may differ from ^^2- 

To elaborate on this discussion, the IR absorption and optical emission PSBs are both 
E A electronic transitions and thus have analogous configuration coordinate diagrams 
(refer to Fig. S3). At low temperatures where only the ground vibronic levels of ^E and ^E 
are populated, the features of the PSBs occur at the harmonic vibronic energies of ^A2 and 
^Ai and are unaffected by the Jahn- Teller anharmonicities of '^E and ^E. Furthermore, for 
typical Jahn- Teller interactions, the contributions of E modes to PSBs oi E A electronic 
transitions are well described by the techniques applied to model Ai mode contributions [7]. 
It follows that a comparison of the PSBs can relate the vibronic structures of ^^2 and ^Ai 
and also the electron-phonon interactions that give rise to the PSBs. 
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In the low-temperature limit, the Huang-Rhys parameter of a transition between the 
ground and n*^ excited electronic states is given by [6, 8] 



where p{u!) is the density of modes of the ground electronic state, Q is the highest mode 
frequency, and (ignoring mode energy shifts) Sn{uj) — a'^{uj)/2huj^, such that a^(a;) is the 
average squared linear interaction parameter of all modes with frequency ou. Likewise, the 
bandshape function due to linear interactions only is [6, 8] 



where Nn is a normalization constant satisfying Nn Jq g{uj) du = 1. Note the minor differ- 
ences in the above expressions to those that appear in Refs. [6] and [8]. 

Using the well-estabhshed expressions of the NV~ electronic states in terms of molecular 
orbitals (MOs) (oi, ex and Cy) [9], the linear and quadratic interaction parameters may be 
derived for each of the NV^ center's electronic levels correct to second-order in electron- 
electron electrostatic interaction (refer to Table S2). Note that only the and ^E' levels 
couple due to the electrostatic interactions and this coupling can be described by the param- 
eter K, [9]. Table S2 clearly demonstrates that the current theoretical model predicts that 
the vibrational parameters of ^Ai do not differ from those of ^A2, whereas the parameters do 
differ at zero-order for ^E and at first- (linear interactions with E modes) and second-order 
(linear interactions with Ai modes and all quadratic interactions) in k for ^E. These results 
are depicted in Fig. S3. As the parameters of ^Ai do not differ from the ground electronic 
level, the above expressions for Sn and gn{i^) corresponding to a transition involving the 
ground electronic level may also apply to a transition involving ^Ai. Thus, the current 
model predicts that Ai and E modes contribute linearly at zero-order to the S and g{uj) of 
the NV" optical emission PSB, but that Ai and E modes contribute linearly at different 
orders in k, (second- and first-orders, respectively) to the IR absorption PSB. 

In summary, the current theoretical model predicts that the mode energies of ^Ai and 
^A2 are the same (i.e. no mode energy shift) and that the contributions of Ai and E modes 
to the optical emission and IR absorption bandshape functions should differ and result in 





gn{uj) = NnS{u)p{u) 



(5) 
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different PSB features. Both of these predictions disagree with the observations made in this 
work. Resolution of this disagreement requires further development of the current molecular 
model of the NV~ center. 



Decomposition and analysis of the vibronic bands 

Adopting the well-estabhshed techniques that were introduced by Maradudin [10] and 
applied extensively by Davies [8] to color centers in diamond, the normalized PSB I{uj) of 
an electronic transition can be described by 

roo 



/oo 
-oo 



oo 

e-^Io{uj) + e-^^—Io^In{uj), (6) 



n=l ^■ 



where S is the Huang-Rhys parameter, 7 is the homogeneous ZPL width, Iq is the homoge- 
neous ZPL shape (as parameterized by 7), g{t) = (l/27r) J^q g{u>)e~'^'^^ duj is the bandshape 
function, Ii{uj) = g{uj) if — Q < a; < Q and = otherwise, In = h® In-i{oj) for n >1, O 
is the highest phonon frequency of diamond, and ® denotes the operation of convolution. 
Note that the individual band components offer the insightful interpretation of being 

representative of all processes of net energy hui involving the creation and/or annihilation 
of n phonons. The total PSB is thus the summation over all n-phonon processes. 

Absorption /absX*^) emission 7em.('^) spectra are related to the normalized PSB via 

4bs.('^0 + <^) hmi^^O - Uj) ji . 

OC — OC /(CU), (7) 

where /icuo is the ZPL energy. The normalization condition -^(1^) (ko = \ ensures that 
knowledge of the proportionality factors is not required to obtain the normalized PSB from 
absorption/emission spectra. Once normalized, the Huang- Rhys parameter can be simply 
evaluated using S — — log(/o), where (/q) is the integrated area of the ZPL. By applying 
Fourier techniques, the following expression for the bandshape function can be derived 

1 /"OO /"OO 

g{^J) = / log{e^e^l*l / [/(w) - h{^J)\e-^^ dw + l}e^'^* dt. (8) 

O J— 00 J —00 

The above expression was found to be non-trivial to evaluate numerically and particularly 
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sensitive to spectral noise. Consequently, it could not be employed to directly obtain the 
optical emission and IR absorption bandshapes. After smoothing the spectra and evaluating 
approximate bandshape functions using the above expression, the evaluated bandshape func- 
tions were then used as seed functions to the iterative process first described by Mostoller 
et al. [11]. The process involves the calculation of each of the normalized PSB components 
via convolutions of a seed Ii{uj) and then using 

oo on 

to obtain an improved estimate of Ii{uj). If a reasonable seed bandshape function was 
obtained directly from (8), we found that the iterative process would converge within only 
a few iterations. 

The fitted IR absorption and optical emission PSBs, together with their individual band 
components /n(^), are depicted in Fig. S4 and Fig. S5, respectively. The plots demonstrate 
that the calculated PSBs reproduce all of the key features of the observed PSBs, but minor 
differences appear towards higher energy. These differences are more pronounced for the IR 
absorption PSB and are likely to be related to the origins of the unexpected phonon energy 
shifts of ^Ai and/or the Jahn- Teller effect in ^E. Further investigation is required. As 
concluded in the article, the features of the optical emission and IR absorption PSBs result 
from electron-phonon interactions with similar phonon modes and that the mode energies 
are greater in ^Ai than in ^^2- It was further concluded that the upward energy shift of 
the phonon modes in the ^Ai level has resulted in a local mode appearing with an energy 
greater than those of the diamond lattice. 

Analysis of the PSBs via a comparison with the phonons of diamond 

Having established the relationships between the PSBs, it is left to explain the origins 
of the PSB features. If one assumes that in ^A2, the NV^ center only slightly perturbs 
the phonon modes of perfect diamond, one can explain the features of the optical g{uj) via 
a comparison with the diamond phonon dispersion curves along high symmetry directions 
([111], [100] and [110]) and the phonon density of states (DOS). As the NV" electronic 
charge is centered on the nearest-neighbor nuclei of the vacancy (NNV), it is expected that 
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FIG. S4. (color online) (a) The normalized observed IR absorption PSB (black) and the calculated 
PSB (red) obtained using the IR absorption §{0:) and (6). (b) The calculated n-phonon components 
(increasing from n=l from left to right) of the IR absorption PSB. The sum of the components 
equals the calculated PSB depicted in (a). 
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FIG. S5. (color online) (a) The normalized observed optical emission PSB (black) and the calcu- 
lated PSB (red) obtained using the optical emission gioj) and (6). (b) The calculated n-phonon 
components (increasing from n=\ from left to right) of the optical emission PSB. The sum of the 
components equals the calculated PSB depicted in (a). 
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the NV center will interact strongly with modes where the NNV undergo the largest relative 
displacement, which occurs when the NNV vibrate directly out of phase. As the NNV occupy 

— * — * — * 

equivalent Bravais lattice sites, this phase condition can be expressed as k- R = tt, where k is 
the phonon wavevector and R is the lattice vector connecting two NNV. For a given phonon 
wavevector direction, the phase condition is satisfied by a single wavevector magnitude. 
The magnitudes corresponding to the [111], [100] and [110] directions are k — and ||, 
respectively, where a is the distance between equivalent Bravais lattice sites in the perfect 
diamond lattice and R is the distance between two NNV. 

As expected. Fig. S6a shows that features of the optical g{uj) approximately correspond 
to features in the DOS and to points along the phonon dispersion curves at which the 
wavevector in the [111] direction has a particular magnitude. Fig. S6 also depicts similar 
comparisons with the [100] and [110] dispersion curves. By considering all of the comparisons 
together, it is clear that one can attribute all of the optical g{u) features to points along the 
phonon dispersion curves whose wavevector magnitudes are related by the parameter i?/a ~ 
1.2. This semi-empirical result is in reasonable agreement with the ratios R/a ~ 1.1 — 1.3 
obtained in previous ab initio calculations [12, 13]. Hence, the elements of the ^Ai and 
vibronic structures observed in the IR absorption and optical emission PSBs can be directly 
related to specific phonon subgroups of diamond. 
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FIG. S6. (a) The comparison of the optical emission g{uj) with the perfect diamond phonon 
DOS p{lo) [14] and dispersion curves of phonons with wavevectors in the [111] direction [14]. (b) 
comparison with the dispersion curves of phonons with wavevectors in the [110] direction (left) and 
[100] direction (right) [14]. Horizontal dashed lines connect related features in adjacent plots and 
the vertical dashed line connects points with the same wavevector magnitude. 
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TABLE S2. The vibrational parameters of the NV~ center's electronic levels as derived using 
the current model of the center. The parameters are expressed in terms of the Coulomb coupling 
parameter k (refer to [9] for further details) and molecular orbital reduced matrix elements of the 
linear and quadratic electron-phonon interactions. The explicit expressions for the reduced matrix 
elements are contained in Table S3. 
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TABLE S3. The reduced matrix elements of the linear and quadratic electron-phonon interactions 

expressed in terms of the NV~ center's molecular orbitals (oi, and Sy). Refer to [9] for further 
details of the molecular orbitals and calculation of reduced matrix elements. 
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